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ABSTRACT 


The objective of this study was to determine by experimental means 
the rotor efficiencies at different radii between hub and tip of a single 
stage axial compressor at its design point, to show the influence of tip 
clearance and end losses. | 

Procedure for calibration and application of pressure probes to survey 
the flow in the compressor were established, and programs were written 
to analyze the measured data. 

Recommendations are made for improvements of the data reduction 
methods, which should precede experiments involving controlled changes 


in the blade tip clearances. 
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I. INTRODUCTION 


The performance of axial flow compressors is greatly affected by the 
blade end and tip clearance losses. No accurate methods exist to predict 
these losses for arbitrary conditions. The main objective of this study 
was to determine by experimental means the rotor blade efficiencies of a 
Single stage compressor at its design point at different radii between hub 
and tip. This efficiency distribution will indicate by how much the so- 
called profile losses are increased by blade end effects. 

At a future date the tests will be repeated with different tip 
clearances to obtain a better understanding of the physical aspects of 
the rotor flows near the blade tips. This study has been undertaken to 
establish the necessary measuring procedure and to set up a data reduction 
method that can be used for these tip clearance tests. i 

The present tests were carried out with a single stage that had 
-solid body rotation bladings designed for forced-vortex flow. It was 
installed in a research compressor of the Turbo-Propulsion Laboratory of 
the Department of Aeronautics at the Naval Postgraduate School. The com- 
pressor has a nominal tip diameter of 36 inches and a hub-to-tip ratio of 
G6. ‘It is driven by a 150 horsepower electric motor at a speed of about 
2300 revolutions per minute, giving a tip speed of 360 feet per second. 

The compressor is throttled ahead of the inlet and discharges into the 
atmosphere through a conical diffusor. The throttling device is installed 
between the inlet bellmouth and the compressor dane. It consists of a 
chamber in which a number of different wire mesh screens or perforated 
plates can be installed to create the desired pressure rise in the com- 


pressor. To operate the compressor along its characteristics, it is 
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necessary to change the arrangement of wire mesh and perforated plates. 
Such changes will be carried out to establish the design point of the 
stage, which is considered to be the operating point where the highest 
efficiency is reached. 

The local rotor efficiencies between hub and tip can be determined 
by taking flow surveys ahead of and after the-rotor at a number of radii. 
Because of the three-dimensional nature of flows in turbomachines, it is 
necessary to determine the velocities and flow angles in axial, peripheral 
and radial directions at each point of a’survey. Such measurements are 
made using 5-hole pressure probes, Type DA 120, of the United Sensor 
Corporation. However, these probes must be calibrated in a flow with 
known velocities and direction at roughly the same Mach numbers as exist 


in the actual compressor. 
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II. EXPERIMENTAL CALIBRATION OF SURVEY PROBES 


A. INTRODUCTION 

The testing apparatus used to calibrate the 5-hole United Sensor 
type DA 120 survey probes is shown in Fig. 1. A Kiel probe was arranged 
as shown to serve as a reference Feeurencnth. Fig. 2 shows a typical 
Kiel probe. Particular pitch angles could be set by use of the apparatus 
shown in Fig. 6. Pitch and yaw angles were defined according to Fig. 5. 
Vertical surveys were restricted to zero pitch angle. 

The mass flow rate from the Allis-Chalmers compressor to the test 


section was controlled by a gate valve. 


B. FLOW FIELD DETERMINATION 

The most accurate probe available is the so-called Prandtl pitot- 
static probe of Fig. 3 which measures true static and total pressure. 
It served as the reference for calibration of the 5-hole probes. Velocity 
profiles were obtained in the horizontal and vertical direction of the 
calibrating duct with this probe. 

For an isentropic process ahead of the probe, the velocity V is, 
from the following any 
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= Cc (tT, - T.) 


For a constant value of the specific heat ratio y , 
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Thermodynamic Process in 
Prandtl Probe 


Introducing the speed of sound a,» at the total temperature qT, » che 
following non-dimensional velocity was obtained: 
x 
Sed 
v 
Pp 
¥-(2(2-(2 a 
iG +t 
i. = P, pa is small compared with the absolute total pressure 


Py , Equation (1) can be expressed by the approximation: 


2s AP) . 
ao(F 4) : 


Surveys were performed in the horizontal and vertical directions in 
the calibrating pipe at various mass flow rates to cover the range of 


velocities that were experienced in the compressor. The readings of the 
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Kiel probe, Pik and Ty , and of the wall static tap in its vicinity, 


Pp » were used as reference values for different mass flow rates and 


W 
inlet conditions. Figs. 9 to 11 indicate a region of uniform velocity 
in the center region of the pipe. A position 4.5 inches from the right- 


hand wall along the horizontal centerline was chosen as the location for 


subsequent 5-hole probe surveys. 


C. CALIBRATION FACTORS OF 5-HOLE PROBES TYPE DA 120 #535 AND #538 


Thermodynamic Process in 
Kiel Probe 


Pat CD (The i T,) 


Similar to Equation (1), with a being the speed of sound for 


tk 


the total temperature obtained from the thermocouple of the Kiel probe, 


aL 
y-1 ? 
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se (ex (es) (3) 
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The mass flow rate at the location of the Kiel probe is obtained from 


V. 


° k 
m=AK,, & v, (4) 


where Kax is a blockage factor, Ok a correction factor for the 


velocity profile if PD, is not the true static pressure and v, is the 


k 
oe : 5 ia a 
specific volume. With Pik ee ae 
seh. FG Tek il 
k B 


—— (5) 
tk (p/P, Y/Y 


Equations (3) and (5) give: 


Vil 3 
a) ee 2/y ia Ne 
Bi vig Fi = fee Gel if (| = $ (6) 
A Kx Se Pix (Yt) Pex Pi 


Since p= Pi, - \, and A was small compared to P Equation (6) 


t tk, ? 


was approximated by, 


(7) 


A similar expression was established for the mass flow rate at the 
location of the Prandtl probe by assuming that the total temperature was 


constant in the pipe, or 


ma JBe/e Tix | ( 4 
ee 2 aL ee (8) 
Se Se Pe 


Since mm and A are equal in Equations (6) and (8) there is 
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Therefore, with the measured quantities P, and P, of the Prandtl 


probe, the total pressure Pik of the Kiel probe and the wall pressure 
DP, near the Kiel probe, it is possible to establish the coefficient 

Y , provided that the above Temaenee data are taken at the same flow 
rate m. The value of y will depend primarily on the Mach number of 


the flow because Reynolds number effects are small in turbulent flows. 


For isentropic processes 


ae 


fw 


where P, » P, are the total and static pressures, respectively, and 

M is the actual Mach number at the pressure ae Hence, for the same 
specific heat ratio y , the Mach number M is only a function of the 
pressure ratio P./D., and is independent of the temperature of the fluid. 


Therefore, it is possible also to consider the ratio 


Ea 7 2 P 
M, = 3 —*-1- = (10) 
tk tk 


as a measure for the actual Mach number M. 

Measurements with the Prandtl and Kiel probes at different flow 
rates, resulting in different values of M. >» gave the relationship 
between. ¥Y and M. that is presented in Fig. 10. 


Rearranging Equation (9) gives the velocity head as 


=e 
12 a 75 (Par Bi) (11) 
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Replacing the Prandtl probe with the 5-hole probe and defining 


P_ - ~p 
iG s 
2 Se (12) 
SE hs 
it follows, 
ie eee 1 
i oe ee (13) 
al i 2 


where PL = P, is the total pressure and Ps is the approximate static 
pressure measured by the 5-hole probe. The location of the pressure ports 
for the 5-hole probes is seen in Figs. 4 and 5.. 


Similar to Equation (10), 


oe (14) 


The relationship between Ky 


15 for the two probes calibrated. These graphs make it possible to 


and M. tSs.SsHown in Bags... J1 12. and 


determine the true velocity head of the 5-hole probe for a measured M 


18 
Different pitch angles were set with the device of Fig. 6. The 

pressure difference between Py and PS is a measure of the pitch 
angle 6. With 

P, -P 

h 

Ky = po > (15) 
ee 


the relationship between this pressure difference and pitch angle is 


shown in Figs. 14 and 17. The relationship between Ka and 6 is 


shown in Figs. 13 and 16. By rotating the 5-hole probes about their 


axes, P, and P,. were balanced and the yaw angle was read directly 


2 3 


from a protractor mounted on the supported frame of the probe. 
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III. ONR THREE STAGE AXIAL FLOW COMPRESSOR 


Figs. 18 and 19 show the compressor installation at the Turbo- 
Propulsion Laboratory. This test rig will be described in detail in a 
report of the Department of Aeronautics, which is being prepared. 

The compressor was designed for a variety of research experiments. 
Six rectangular ports are arranged in the upper casing to accommodate 
the survey probe carriage shown in Fig. 20. The location of these ports 
makes it possible to survey circumferentially in a iiss arc at positions 
ahead of and after the rotors and stators. 

The rotor and stator blades are removable and during this study only 
the first rotor and third stator were installed. The relative position 
of the rotor and stator for this "expanded" single stage configuration 
ts seen, in Figs 2l. 

The stator and rotor blade profiles of the solid body rotation stage 
are Saiiiyini Fig. 22. The blade chords varied from 2.00 to 3.40 inches. 
The solidity of the rotor varied from 1.150 to 0.903 from hub to tip and 
that of the stator from 1.039 to 0.850. 

Inlet guide vanes (IGV) are located ahead of the rotor and two 
rows oe exit guide vanes (EGV) are located after the compressor stage. 

There are thirty rotor blades per row and thirty-two IGV's, stator 
blades, and EGV's per row. 

A more detailed description of the blading is given in Ref. 2. 

As shown in Fig. 23, Kiel probes were located in the survey holes 
along the length of the compressor, which measure the total pressure at 
fixed radial positions. 

The inlet bellmouth was surrounded by a screen enclosure to reduce 


disturbances by wind gusts. 
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IV. INSTRUMENTATION 


Wet and dry bulb mercury glass thermometers were located on the 
inlet screen to determine the relative humidity of the air entering the 


for different 


bellmouth. A program for determining the gas constant Ra 
humidities is given in Appendix B. | 

The pressure drop across the inlet nozzle was measured with a Meriam 
‘Micromanometer, seen in Fig. 24, which had a reading accuracy of 0.001 
inch of water. 

Kiel probes were located in the survey holes and static pressure 
taps were located in the same meridional plane. These pressures were 
measured on the large water manometer board of Fig. 25. 

A United Sensor 5-hole probe was located in the survey carriage, 
seen in Fig. 20, and the various pressures were read from the U-tube 
manometer board of Fig. 25. 

Forward and aft bearing temperatures were read by a Brown Electronik 
temperature indicator, as seen in Fig. 26. 

A BLH SR-4 Torque Pickup type A, with a maximum torque of 5000 in.-1lbf. 
and a maximum speed of 5000 RPM, was attached directly to the compressor 
shaft as shown in Fig. 27. The torque was measured by a Doric Integrating 
Microvoltmeter, Model DS-100, after the signal had passed through a 
d.c. amplifier, and a Krohn-Hite Model 3750 Filter, which are shown in 
Fig. 26. The calibration of the torque pickup is described in Appendix 
A. 

The rotational speed was measured by a flux cutter located at the 


shaft end. The readout was displayed on the Dynasciences Electronic 


Counter of Fig. 28. 
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The barometric pressure was recorded with a Princo Fortin Mercury 
instrument. Gravity corrections were made according to the instruction 


booklet and temperature corrections were made according to Ref. 3. 
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V. DESIGN CONDITIONS 


A. DERIVATION OF PERFORMANCE PARAMETERS 

To determine the design point for the ONR Compressor it was necessary 
to establish parameters that give a measure of its performance. The 
performance was based on the conditions that see at the mean radius 
of the compressor, 14.4 inches from the centerline. As stated in Ref. 2, 
the solid body blading has fifty percent reaction at the mean radius. 
Figs. 29 and 30 show velocity diagrams for such bladings. 

The efficiency was defined as the ratio of the theoretical power 


P. and the actual drive power P, , or 


th d 
ie 
th 
= (16) 
Fa 


The thermodynamic process through the compressor stage is represented 


in the figure below. 
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The inlet conditions are those ahead of the rotor and discharge conditions 
are those after the stator (see Fig. 30). 


The theoretical power is 


P,, =m ME. (17) 


From the H-S diagram, assuming y is constant, 


a ! a3 
Py =n c. (Ti T,5)J (18) 
where Teg! is the discharge total temperature for an isentropic conm- 
pression, and Ts is the total inlet temperature. For isentropic 
conditions and AP. = Pig - Piss small compared with the absolute total 
inlet pressure Pi > there is 
AP 
te -1 aoe 
Wear a ta y AS eo) 
ti ; 
AP. R 
a pe pe } = _G& ps asd 
where F=1 By oe With C F vo? 


Equation 18 becomes 


AP, : 
Pin m Ra Tey eis (20) 


The mass flow rate was determined from the pressure drop in the inlet 


pellmouth Pry = P, » Where P is the atmospheric pressure and 


a” Py ta 


Py is the wall static pressure at the minimum diameter of the bellmouth. 


With Bernoulli's equation, Ref. 1, page 203, 


26 Bei T 
Vy Pe ae (op,) (21) 
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The mass flow rate is 
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where Ke is a blockage factor which is taken as wnity for this derivation. 


With Equations (21) and (22), Equation (20) is 


Jk AP 
estre2 tA t 
Fn “Ey J728o pO Wy Fe Tea BS (23) 
tA ti : 
The actual drive power Ps is given by 
-yMum 
P= M 30 (2h) 


where M is the torque and N the rotational speed in RPM. 


The flow coefficient g is defined according to Ref. 13 


(25) 


where v. is the average axial velocity and U the peripheral speed at 
the mean radius. 


The work coefficient is defined as 


Fpaeee (26) 


U/e 


For a constant mean radius there is from Euler's turbine equation 


(Ref. 1, page 259). 
T= (27) 


where wv, is the change of the peripheral components of the absolute 
velocity in the rotor. 


The so-called head coefficient is defined as 


AH. 
y = oes (28) 
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From Equation (25), (26) and (27) 


B. EXPERIMENTAL PROCEDURE 

As described in Sections III and IV, Kiel probes were used in deter- 
mining the total pressures. The static pressures were taken as the wall 
pressures. From Fig. 30, there is V. = V. at the mean radius. Hence 
P43 - Pi = P2 - Py where Pie and Pia are the total pressures 
measured by the Kiel probe and Py and P2 the wall static pressures. 
Therefore the performance of the stage could be based on either total 


or static pressure variation. The static ‘pressure approach was used. 


¢ 


C. DATA REDUCTION PROGRAM 

Appendix C shows a card See aus the Monroe Calculator Model 1656, 
which determines the performance parameters from measured input data. 
Numerical results were obtained for twenty-one throttling configurations, 
with three different readings taken for each one. These data are listed 


in Table C-I. 


D. RESULTS 

The efficiency, and the work and head coefficients, are shown in 
Bigs 7.\31,. 32 and 33 as functions of the flow coefficient ». fhe 
linear interdependence of +t with g in Fig. 31 is characteristic for 
most axial compressors, as os in Ref. 1, page 352. Near the surge 
point of the compressor the slope of the curve changes. Two configurations, 
wire mesh #4 and perforated plate A and wire mesh #5 and plate A, pro- 
duced conditions close to surge. The efficiency and head coefficients 


also showed substantial changes for these configurations. 
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A least-squares subroutine, LSQPL2, from the IBM 360 System Library 
was used to determine the best slope through the data of Fig. 31. 

The maximum efficiency was about 92.5 percent. The throttling 
configuration which produced this highest efficiency was that of wire 
mesh #3 and perforated plate A. Figs. 31 and 33 show that this con- 
figuration operated sufficiently far from the surge point. Hence the 


flow surveys were carried out for this throttling arrangement. 
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VI. ROTOR SURVEYS 


A.- SURVEY METHOD 
The United Sensor 5-hole probe Type DA 120 #535 was used in the 
survey carriage of Fig. 20. Surveys could be made over an are of 15° 


in circumferential direction and from hub to tip in radial direction. 


B. DATA REDUCTION METHOD 
1. Corrections for Ambient Pressure and Speed of Rotation 
The surveys were made over an extended period of time at dif- 
ferent atmospheric conditions and rotor speeds. Changes in the ambient 
temperature Tho » and ambient pressure Pia > and compressor speed N 


have an effect on the velocities and pressures at the survey stations. 


These changes are reflected in the fixed reference total pressure, 


Pur? 


and Pp the reference static pressure, measured at the locations shown 


in Fig. 34. ‘The quantity dp = Pip - Pp is a measure of the velocity 
head of the flow through the compressor, even though Pip and Pp 


were not measured at the same station. 

From the test data, sample calculations shown in part 1 of 
Appendix D indicated that an average density change of four percent 
occurred in the machine. Thermodynamic relations were therefore used 
to calculate the velocities and densities at the survey stations. 
However, to apply ponnee eiense aor changing values of P ; Tea and 
N it is sufficiently accurate to assume incompressible conditions. 

1s was taken as the total pressure ahead of the compressor 


tI 


blading at a particular radius, and Pip as the total pressure down- 


stream of the compressor blading. For a fixed throttling configuration, 
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the pressure drop from the ambient pressure P to P can be 


expressed by 
= £y- 
By, te beet SC V (29) 


where i is the axial velocity through the blading. From the compressor 
discharge to the diffusor exit, where the pressure Pin exists, the 


pressure change is 
= Lye 
Ba.) iP c V (30) 


Cy and Ch are constants for a particular blading and inlet nozzle at 


a fixed throttling configuration. Combining Equations (29) and (30), 


P = Ly? 
P Pir (Cc, +c Vv (31) 


tD 


From Euler's turbine equation for incompressible flow 


bay = ee = uy (32) 


V 
where Y is the pressure rise or head coefficient. With = = 3 


and Equations (31) and (32) 
, a i) 2 


It can be eer that a compressor with fixed geometry operates at the 
same operating point independent of P = Thy and Ui. 

Similar to Equation (29) and assuming a perfect gas, Pi, - Pip 
can be expressed by 


P 
tA a 
P - P = (¢ See ys 
tA “tR OR 2 RT 6 


where Ce and o are constants. With U = fal 3 
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where Kp is a constant that is not dimensionless. Since the wall 


reference pressure Pp is lower than Pip 5 
ae iy 
Puno R = 2a 


since the loss in total pressure from the Kiel probe to the wall tap 
and the velocity head are multiples of £y° ‘ 
With the equation for o , 
12 


tA ae 
as 
c Tta & 


or, because @ and € are constants, 


D ‘ 

N 

= - — (34) 
PLA R 5 ER, = 

A similar approach can be carried out for the total pressure Pia ; 

and the static pressure Pi >» measured by the flow probe at station l. 


Hence similar to Equations (33) and (34) 


Bg Seat | N = es 

ee fe be Se) 
Pha = VER =) 

a es 


al! N 
a 


ash Pip and P, are measured at the same speed N, temperature 


T., » and pressure P,, , there is from Equations (33) and (35) 
Je = P. K 
— —- = — = K = constent (37) 
ta ftp *R 
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The purpose of this approach was to correct the measured pressure 
Pia to a constant average ambient pressure (P,) - With this average 
O 


pressure corresponding to an average total reference pressure (Pp) 5 
O 


an average speed No and an average temperature (q,) » the so-called 


corrected pressure (P,) that exists at (P,) a1, 4 (t, ) is obtained. 
c O 
From Equation (37) 


(2B. = (PL?) 

~ Wipe Sa SST 
Pp = (P 7 Paes 
moe tR’y tA tR 


and rearranging, 
Gap 4 ) el 
P = (2 =- (P. = P (38) 
tl a A 0 A tl Pia - Pip 


In asimilar manner, using Equations (34) and (36) and introducing 


= 


the average static pressure (p,) » there is 


(Pip) 2 (835) 


e ae > 


Since (P, is obtained from Equation (38) the corrected static pressure 


1) 
e 
aGescacaon 21s kniown: 


Equations (38) and (39) can alsobe used to establish the corrected 
total and static pressures at station 2 after the rotor, by replacing 


index 1 by index 2. 


2. Performance of the Rotor 
With these corrected pressures it is possible to determine the 
performance of the rotor. 


To obtain the velocity Vy at station 1, the total temperature 


a was taken as the average ambient temperature (T, ) . Then by Fig. 36, 
O 


tl 
. a aye 
with Ch = Ra g Wied. 
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where * is a dimensionless quantity. Pia and P, are the corrected 


pressures, and (R,) is the average gas constant obtained from the 
0) 


average total temperature. The density Py is obtained from 
nce 
eR 
and (By) 
O O 
Then with Ta = (T,) ’ 
L 


v (R.) 
Re a, (=) Fel de YO (42) 
= > ae Oe / 
t Cy) Pel Fy) nie 
The total temperature T after the rotor is determined from the 


t2 


average torque My » the average speed No and the average flow rate 


w . Since Tho - Thy = Tio 


calculations of Appendix D, it is sufficiently accurate to determine the 


- (T, ) was small, as found in the sample 
(0) 


flow rate w by 


(42) 
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where [Py is the average pressure drop in the inlet bellmouth at its 


smallest diameter Dy é 


Then with My in ft-lbf and No in RPM, the actual driving 


power is 
ie 
Fimiueg “0 30° * Cp ‘fia %)) 
T 
+2 
= w(R,) CEs Shel 
Ge T), 
T TiN /30 
or ey =i + ee (43) 
Arp Wg) yer Ow 
0 
Then, similar to Equation (40), 
y-1 
2 Bay Yo iiee Ra 
Payee ogee ea . Mb 
Vo TaN ie (r=) Ene (an Ch.) 
) 0 
and by Equation (41) 1 
Po Me) Ce) aa (45) 
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Fig. 37 shows the measured pitch angle \A‘' and the measured yaw 
angle a’ . The peripheral velocity a > axial velocity Mee and 


radial velocity ee are, 


es = V*% cos X' cos a’ 


Ve sin -x' 
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Vw = Vaecos: sin a" 


31 


Ue f m | eT o> 
WD) ABU SS Sh MIA ce Sy 


ios . oi 1 
Lat 7 Gyr? q” sie EF re ovicd 


Wie aft PP 


— a ee oe a 
La 4 
ore *, 
» + y : a7 
wy wn © *; 


> +) 


-(08) aotiacp® of teLimta ose 


t a ry f T ' Bs ' . —_ 5 — 
a -~ ° faded ¥ 
x , =e ; ) te —_ - | Ti ” rv 
aJ \ allt J \s >> ~"7 -_ 
. q ie - 
P (iH) wottaupt ed f 
ly Se) ( 2) \oa> VP 
r + = | ck 
aptneiuiaiimanes seal w 
ie oe Ci \cet lest? oe 
Ww : d “A : avy SJ ae 


7 

Se r<4 it 6 '< sfeome ctotte beruveasem af? ewode YE saht = 
é ‘cf? “2 . re 

: y stooley Ieixe ,. * Pr | yilicaolev jatsdqPreq emf . “® otpme 


,o% ag? _Wtooler 3 
elem ¢ 
: : “ay he 
‘a eon “f sos *Y = #2 roe ~—_ 
S “a la 
yo i at 


i nle Wie Oy . a 
oi , ie 4 , 7 


. ‘4 7 
'o dia ‘Ss sop *Y = *7 


and the relative velocity w* is 


2 2 ai 
W* = ¥ L *% ¥* * Te 
[(U v. yr + Nerve ] 


where Ux g 
g Ra Ty 
0) @) 
The actual flow angles are 
Vv * 
: U ? 
sin a = =z = sin @' cos i' 
U*-V * 
3 = u 
Sony 6 = art 
Vv * 4 ; 
eine 7 _ Sin r 


V* cos a cos @ 


Velocities denoted with an asterisk are corrected to average ambient 
conditions and average rotor speed. 
Because the flow properties vary in the peripheral direction, as 


illustrated in Fig. 39, it is necessary to obtain average values by 


integration. From the moment of momentum equation 


M = rb - f dm, R, V (46) 
where M, is the axial moment. With R, = R, = R 
M J. Me Vo J. ah Vy 


ay Ch /rmcemmre 

ie cary pal a ee ee (47) 
(2) (1) 

The mass flow at a particular radius R for a differential dR, is 


dh =p V RdedR 
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where 7m are the mass averaged peripheral velocities. 
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The power necessary to drive the blading is 


P= AiC, Of, = dh (R,) ea (50) 
Rr 
rotor en | SS 
a(R T e(R T 
EO GEG on 6 
Equation (50) becomes 
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The mass averaged total pressures at particular radii, with the 


corrected dimensionless velocities, are 


* * 
FE Py a ey O° 


£,, = (52) 
ue J p,* La dé 
(1) 
and 
mw ps = Viz PL ae 
Pre = i = “4 a (53) 
(2) 2 e2 
With the values of Biss and a thus determined, Fig. 36 gives 
Rapes fs 
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This relation was used to calculate the rotor efficiency at particular 
radii. 
C. COMPUTER PROGRAM 

A sempatee program was set up in Fortran IV to determine the 
rotor efficiency with the IBM 360 Computer. 

Appendix D lists the aeoeea: with computer printouts for various 
radial positions. 
D. RESULTS 

The main result of this study is shown in Fig. 38, which is a 


plot of the calculated rotor efficiency as a function of the radial 
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position. At the mean radius the efficiency is about 94 percent. Fig. 
38 shows an increase of efficiency toward hub and tip. From these peaks, 
the efficiency drops radically toward the outer and inner walls of the 
channel. 

The wakes of the inlet guide vanes were quite prominent ahead of the 
rotor. As seen in Fig. 39, they were not detected behind the rotor, 
since the pressure probe gives a time average of the periodic flow. 

The pitch angle of the flow, and thus the radial component of the 
velocity, was small near the hub, but increased steadily toward the tip. 
At a radius of about 17.5 inches, however, the radial flow reversed 
itself. 

Near the hub region of the blade, the axial velocity increased 
across the rotor, whereas it decreased near the tip region of the blade. 
This agrees with the test results of Ref. 2. 

The peripheral component of velocity, Vu , increased with increasing 
radius up to about 17.5 inches. This was expected since the blades 


were designed for forced vortex flow. 
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VII. CONCLUSIONS AND RECOMMENDATIONS 


The radicaldrop of the efficiency near the hub and tip of Fig. 38 is 
caused by blade end and tip clearance effects. However, it is very 
unlikely that the actual rotor efficiency can exceed unity. It is 
necessary, therefore, to evaluate the effect of the assumptions that are 
basic to this study. 

For the determination of the rotor efficiency it was assumed that the 
flow had cylindrical stream surfaces. The fact that the pitch angles 
increased towards the outer casing makes this assumption inaccurate. 

The surveys were made with the 5-hole probe. Calibration tests were 
made to determine the true static pressure and the pitch angles. Through- 
out this testing it was assumed that Pia was the true total pressure. 

At pitch angles of 20° this may not be true. The stem effects were also 
ignored. 

It is recommended that further calibration tests be performed with 
the 5-hole probe to determine the effect of pitch angle on total pressure. 
Since pitch angles exceeded 25° in the compressor, it is necessary to 
cover a wider range of pitch angles in these tests. 


To establish the magnitude of the tip clearance and end losses it 


is necessary to perform tests with increased tip clearances. 
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APPENDIX A (Torque Calibration) 


The calibration of the torque readout instrumentation, which is showm 
in Figs. 26 and 27, was carried out by applying a static torque of known 
magnitude. This procedure consisted in hanging different weights from 
a twenty inch moment arm that was attached to the drive shaft, and 
recording the counts of the readout instrument. Table A-1 is an example 
of one of the three torque calibrations made. The data from these three 
calibrations is shown below. This curve produced a calibration constant 


of 17.3333 in. lbf per count. 


3000 
Torque 
(in. 1bf) 


2000 


1000 


25 50 75 2 160 125 


Instrument Counts 
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TABLE A-1 


EXAMPLE OF TORQUE METER CALIBRATION 


METER READING 


IN COUNTS 
TOTAL WEIGHT APPLIED TORQUE 


IBF IN-LBF 
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APPENDIX B (Gas Constant Variation) 


The gas constant R is defined as 


G 


_ 1545.43 |ft-lbf 


Ra M Tea=on _  (B-1) 


where M is the molecular weight of the atmosphere which is a mixture 


of air and water. 


For gas mixtures, the number of moles 


a Amount of Constituents in lbs 
ro M, of Constituents 


The molecular weight of a mixture of air and water is 


at n, (28.97) + n (18) 


W We Gee 


where n, is the number of moles of dry air, with M = 28.97 and n 


h 
is the number of moles of water vapor, with M=18. 
ise 1545.43 (n, + n) BB 
sne(20.o7 em) ; § 
G n, (20.97) +n, (1 
Using the wet and dry bulb temperature measurements, ee and 


Te , obtained at the compressor inlet, the relative humidity and the 
vapor content in grains per pound of dry air for 100 percent relative 
humidity is determined from a psychometric chart [Ref. 3]. With 7000 


grains per pound, 


_ (r) @grains/1b dry air ) 
> 7000 


where r is the relative humidity. Equation (B-2) is solved with the 


card program of Table B-1 for Monroe Calculator Model 1656. 
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APPENDIX C (Mean Radius Compressor Performance Program) 


The performance parameters of Section VI are calculated on a Monroe 


1656 calculator with the program of Table C-3. 


The input parameters are: 


~ RPM 


Atmospheric Pressure (in. Hg.) 
Pressure drop across enclosure to inlet (in. H,0) 


Total inlet temperature (°R) 


Gas Constant =) 


Pressure drop across inlet nozzle (in. H,0) 

Pressure rise across stage (in. H,0) 

Pressure difference just ahead of the rotor (in. H,0) 
Torque counts 


Compressor Speed (RPM) 


With these inputs, the parameters $, 71, 1) and ¥ are calculated. 


To find the design point, 21 throttling configurations were tested. The 


results are listed in Table C-l1. The dimensions of the wire mesh screens and 


perforated plates are listed in Table C-2. 
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TABLE C-2 


THROTTLING MEMBRANES 


(1) SCREENS 
1 2 
Mesh Number (mrt) y 10 
Wire Dia. (<5 IN) 28 iG 
(2) PERFORATED PLATES 
A B 
Hole Dia. (IN) 3/4 5/8 
Spacing (IN) at 7/8 
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APPENDIX D (Data Reduction Program for Rotor Surveys ) 


A. SAMPLE CALCULATIONS 
1. Density Changes 
Reference 4 shows, that for one stage operation, the pressure 
fise is about je in. H,0, for a volume flow rate VW of 1000 ft>/sec. 


and a drive power of 135 H.P. At P,, = 14.7 psia and T,, = 520R, 


tA tA 


the weight flow w is 


w= Wop 


76.3 lbm/sec 


Then the maximum possible temperature rise of the air is, with 


a Btu o 
es = wok nee 


The average axial velocity at station 1 of Fig. (34) is 


Vi= ee 221 ft/sec 
a A 
For flow angle ao = 25° . 
Na 
= = 2 / 
vy SiGe a 3 ft/sec 
For an isentropic expansion from Pip ; ThA to Py > qT) 
‘ > 
Sed. = ‘Ib en 
2gJC ta al 
iE 
Al 
or — = 0.9906 
tA 
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Pp 
| 73 = 0.9673 
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The survey tests were carried out at an average P 


a of mh psi. 


Hence p p ip 

L 1 ta 

Fr See pe he 
a tA us 


After the rotor the average total pressure Pio is 14.84 psia. Since 


e ‘ ae 
Viz = Vap and assuming a, = hor, 


p22 
2 cos QL, 


= 288 ft/sec 


: T 
a 252 020669 


0.9548 


therefore LD pe 0.9%70 
=> 
B. DATA REDUCTION PROGRAM 
The data reduction method developed in Section VI is performed 
by the IBM 360 Computer. Measured pressures and temperatures are cor- 
rected for average conditions. The program is listed in Table D-l. 
The velocities ahead of and after the rotor, together with the rotor 
efficiency for several radial positions, are listed in Tables D-2 to D-15. 
i. Data Imput 


Data are read in as follows, first for the condition ahead of 


and then for those after the rotor. 
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FORTRAN NAME 


AHEAD 


OF ROTOR 


DELPAE 


PERL 


AFTER 
ROTOR 


Ne 


PAHG2 


1 


DLPIR2 
DLPAR2 
DEPE2 
DLPN2 


DLK3A2 


DLKTA2 


DLK5A2 
DLS 3A2 


DLPAE2 
RG2 


Ere 


PER2 


FORMAT 


Lye 


FLO.5 


FLO25 


BLOL5 
FLO.5 
BLO. 
PLOND 


F10.5 


F10.5 


F10.5 


BLO 


F10.5 
F10.5 
F10.5 
F10.5 
F10.5 
F10.5 
F10.5 
F10.5 
F10.5 


F10.5 
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DESCRIPTION 


Number of Positions Surveyed 
Circumferential Direction 


Atmospheric Pressure (in. Hg 


in 


-) 


Temperature used in correcting 


atmospheric pressure ( F) 


Pressure difference P.47Pp (in. H,0) 


Pia "Pe (in. H,0) 
Pay -Py (ing H,0) 
Py (in. H,0) 


Pia - Kiel total pressure ah 
rotor (in. H,0) 


Pe = Kiel total pressure ah 
tA i 
IGv (= Pip) (in. H,0) 


ead of 


ead of 


rad 


Pia > Kiel total pressure after rotor (in. H 
P,, - static wall pressure at rotor 


(= PQ) (in. 4,0) 


PLA ~ enelleetrne (in. H,0) 
EASE 
1bm-9R 


Ambient Total Temperature, T 


Gas Constant ( 


fe) 
d@ , used in integration ( ) 
Yaw Angle a (°) 
Rotor Speed (RPM) 


Torque (counts) 


ie PL-Po 
fe ee 


x 
P= P. (in. 1,0) 


Circumferential position @ 


eer 


eo) 


0) 


Listing of Rotor Performance Program 
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Figure 4. Enlargement of the Tip of a 1/8" Diameter 
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Figure 17. K, Versus Pitch Angle for 5-Hole Probe DA 120 #535 
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Figure 20. Survey Carriage 
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Figure 24. Meriam Micromanometer 
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Figure 33. Head Coefficient Y¥ Versus Flow Coefficient 9 
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Figure 35. Blading at a Particular Stream Surface 
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Figure 36. Temperature - Entropy Diagram of the 
Thermodynamic Process Across Rotor 
Blade (B = 2gJC,,) 
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Figure 37. Velocity Diagram Showing Relationship 
Between Measured and Actual Flow Angles 
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Figure 39. 


Velocity Before and After Rotor Versus 
Circumferential Position at R = 15.0 in. 
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